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Summary. — We develop a new method for describing the dynamics of 3- 
dimensional thermal plasmas. Using a piecewise constant 1-particle distribution, 
we reduce the Vlasov equation to a generalized Lorentz force equation for a family 
of vector fields encoding the discontinuity. By applying this equation to longitu- 
dinal electrostatic plasma oscillations, and coupling it to Maxwell's equations, we 
obtain a limit on the magnitude of the electric field in relativistic thermal plasma 
oscillations. We derive an upper bound on the limit and discuss its applicability in 
a background magnetic field. 

PACS 52.65.Ff - Fokker-Planck and Vlasov equation. 
PACS 52 . 35 . Fp - Electrostatic waves and oscillations. 



1. Introduction 

High-power lasers and plasmas may be used to accelerate electrons by electric fields 
that are orders of magnitude greater than those achievable using conventional meth- 
ods [1]. An intense laser pulse is used to drive a wave in a plasma and, for sufficiently 
large fields, non-linearities lead to collapse of the wave structure ('wave-breaking') due 
to sufficiently large numbers of electrons becoming trapped in the wave. 

Hydrodynamic investigations of wave-breaking were ffist undertaken for cold plas- 
mas [2, 3] and thermal effects were later included in non-relativistic [4] and relativistic 
contexts [5-7] (see [8] for a discussion of the numerous approaches). However, it is clear 
that the value of the electric field at which the wave breaks (the electric field's 'wave- 
breaking limit') is highly sensitive to the details of the hydrodynamic model. 

Plasmas dominated by collisions are described by a pressure tensor that does not 
deviate far from isotropy, whereas an intense and ultrashort laser pulse propagating 
through an underdense plasma will drive the plasma anisotropically over typical acceler- 
ation timescales. Thus, it is important to accommodate 3-dimensionality and allow for 
anisotropy when investigating wave-breaking limits. 

Our aim is to uncover the relationship between wave-breaking limits and the shape 
of the 1-particle distribution /. In general, the detailed structure of / cannot be re- 
constructed from a few low-order moments so wc adopt a different approach based on a 
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particular class of piecewise constant 1-particle distributions. Our approach may be con- 
sidered as a multi-dimensional generalization of the 1-dimensional relativistic 'waterbag' 
model employed in [5] (for a discussion of the relationship between our approach and [5] 
see [9]). 

We employ the Einstein summation convention throughout and units are used in 
which the speed of light c = 1 and the permittivity of the vacuum eg = 1- Lowercase 
Latin indices a, b, c run over 0, 1, 2, 3. 

2. Vlasov-Maxwell system 

Our attention is focussed on plasmas evolving over timescales during which the 'dis- 
crete' nature (collisions) of the plasma electrons can be neglected and the plasma ions can 
be prescribed as a background. Such configurations are well described by the covariant 
Vlasov-Maxwell system [10, 11] which, for the purposes of this paper, is most usefully 
expressed in the language of exterior calculus (see, for example, [12,13]). We will now 
briefly summarize the particular formulation of the Vlasov-Maxwell system employed 
here. 

Let {Ai,g) be a spacetime with signature (— , -I-) for the metric tensor g. Each 
point p e is associated with a space £p C TpJ\4 of future-directed unit normalized 
vectors on Ai, 

(1) £p = {ix{p),x) e TpM : gab{x{p))x''x^ = -1 and x° > 0}, 

where gab arc the components of the metric 5 in a coordinate system (a;") whose patch 
contains p and (x°, x^) are induced coordinates on TM. The total space £ of the bundle 
{£,Ii,M) is the union of £p over p <E and 11 is the restriction to £ of the canonical 
projection on TAi. 

Naturally induced tensors on TM include the dilation vector field X 

(2) X = x'^dY, 

the vertical lift ★l^ of the volume 4- form +1 from A4 to TA4 and the horizontal 4- form 
#1 

(3) #1 = v/|detfl|^da;°^ A dx^" A dx^" A dx^" 

where dx"'^ is the horizontal lift of c?x" from M to TAi (see Appendix A for further 
details) and q = (gab) is the matrix of components of g. 

The Vlasov-Maxwell system for / (a scalar field on TA4 whose restriction to £ is the 
plasma electron 1-particle distribution) and the electromagnetic field F may be written 

(4) Lf ~ 0, 

(5) rfi^ = 0, d*F = q*{Nion- N) 

where ~ indicates equality on restriction by puUback from TM to £ and 

(6) L^x'^id^ + C)eTTM, 

(7) fa - --F'adb e TM 
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with TO the mass and q the charge of the electron [q < 0) and F";, = g"''^Fcb the com- 
ponents of the electromagnetic 2-form F = ^Fabdx"' A dx^ on M. The metric dual of a 

vector V is defined by V{Y) = g{Y, V) for all vectors Y on Ai and ★ is the Hodge map 
induced from the volume 4-form +1 



(8) ^1 = Vidctfll dx" A dx^ A dx'^ A dx^ 

on A4. The components of the electron number 4-current N = N°-{x{p))da at p € A4 
are given as an integral over the fibre £p = Il~^{p) 



'n-i(p) 



(9) 



n-i(p) goc(a;(p))a;^ 



and the ion number 4-current TVion is prescribed as data. 

The measure on £p in (9) is induced from the 3-form lx^I, 



-V 1 



#1 = ^d^ ^x'^eabcddx"" A dx^" A dx''" 



3! 



(10) c^-^}^^dx'" Adx^" Adx'" 



9oc 



where Cabcd is the alternating symbol with eoi23 = 1- 

The Vlasov-Maxwell equations constitute a non-linear integro-diffcrcntial system. Di- 
rect calculation of its solutions for general plasma configurations is difficult and, to pro- 
ceed analytically, it is common to approximate the above as a finite number of moments 
of / in satisfying a non-hnear field system on M (a so-called 'fluid' model). How- 
ever, there are difficult issues associated with closing the resulting field system (see, for 
example, [14]) so we opt for a different approach. Our strategy is to reduce the system 
by employing a discontinuous /, and to proceed we need to cast (4) as an integral. 

One may rewrite (4) as 

(11) d{fu)^0, 

(12) Lo^LLi^l^ Aix#l) eAeTM. 

Integrating (11) over a 7-chain A C £ and applying Stokes's theorem yields 



(13) / fu^O, 

J OA 



with dA the boundary of A. For differentiable distributions, this equation is equivalent 
to (11); however, since it makes no reference to the differentiability of /, it may be 
regarded as a generalisation of (11) applicable to discontinuous /. 
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3. Evolution of discontinuities 

Equation (13) may be used to develop an equation of motion for a discontinuity, which 
we choose as a local hypersurface 7i. Suppose that A in (13) is a 7-dimensional 'pill-box' 
straddling Ti. We may write dA = <t+ + (T_ + (Tq where (t+ and <t_ are the 'top' and 
'bottom' of the pill-box and ctq is the 'sides' of the pill-box. Thus, in the limit as the 
volume of A tends to zero with ct+ tending to a and (t_ tending to —a, we recover the 
condition 

(14) [f]a*Lo = 0, 

where the image of tr is in and [/] = a\j -I- <7*_j . Thus it follows that a finite 
discontinuity in / can occur only across the image of a chain E satisfying 

(15) ITu = 0. 
Suppose that E may be written locally 

E : V X P ^ £• C TM 

(16) (x^$^C')^(a;^i'' = E''(x,e)) 

for V C M., where E'' denotes the component of E, and E T) C M^. It is then 

possible to translate (15) into a field equation for a family of vector fields on V given 
as 

(17) V^ip) = V^{x{p))da = t'^{x{p),i)da 

where, since gab{x{p))x°-x^ = —1 at p E £ , it follows 

(18) g{V^,V^)^-l. 
Using (12, 15) it follows 

(19) E*(ti*l'^Atx#l+*l'^AtL/,A-#l) = 0. 

V ' V ' 

(a) (b) 

Consider first the term (a) in equation (19): 

(20) E*(ti*l^ Atx#l) = E*(i"tg« *1^ A tx#l) 

where (6) and ifv * 1^ = have been used (see (A. 3) in Appendix A). Thus 

(21) E*(i''ia« * 1^ A tx#l) = *Ve A E*iA'#l, 

since lqh * 1^ = {gat * dx^)^ (see (A. 10) in Appendix A). Furthermore, using (A. 5) in 
Appendix A, it follows 

(22) E*(dx"-^) = DV^ + dt" 
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where D is the exterior covariant derivative on A4 and d is the exterior derivative on V. 
Using (10, 22) it foUows 

(23) = V\d^\^eabcdV^"{DVl + dt") A {DV^' + dt^) A {DV( + dt") 
and (20) is 

(24) E*(ti ★ l'^ A = A ^^^^\A^\eabcdV^DV^ A dS= A dt'' 

since ^Vj A DV^ A L'T^'' = (dini(Al) = 4) and dS" A dtfi /\ dt" = (dim(2?)=2). 

Since -kV^ A DV^ -(Vy^F^)'' ★ 1 and \/\detg\eabcd = ^9, * 1- it foUows (24) 

may be written 

(25) E*(iL ★ 1^ A Lx#l) = %A Vv,% A A d^^ A 
where the family of 2-forms Cl^ on V is 

(26) = -g^-Q^gacgbddx" A dx'^. 

The second term (6) in (19) can be rewritten using a similar procedure : 

(27) E*(*l'^ A iLixm) = *1 A E*(/,iix#l) 
and from (6, 7, 10) it follows 

(28) iLOx#l = -^-V\d^^F''Yi''x'eabcddx^'' A dx'^'f'. 

2! m 

where (A. 9) and (A. 10) have been used. Then 

(29) E*(tiix#l) = -l^^ld^lF'^V^-VleabcdiDV^' + dt'') A (W/ + dE'^) 
and 

(30) E*(7^1^ A LLix^l) = - — % ALv^FAflfA A d^^. 

TO 

Combining (25, 30) and (19) yields 

(31) E*w = i/^ A (Vv^Fe - ^ty^F) A A d^^ A df = 0. 

Acting on (31) successively with iy^, ig/g^i and Lg/g^2, and noting that ty^ri^ = and 
ff(V^?,F?) = -l, yields 

(32) (Vy.V^f-— ty,i^)Al7. = 0. 

' TO ' 
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Thus, solutions to (32) may be obtained by demanding that is driven by the Lorentz 
force 

(33) Vv,% = -LV,F. 

However, ahhough (33) is simpler than (32), there are simple solutions to (32) that do 
not satisfy (33); we will return to this point shortly. 

4. Non-linear electrostatic oscillations 

A laser pulse travelling through a plasma can excite plasma oscillations, which induce 
very high longitudinal electric fields. Due to nonlinear effects, there is a maximum 
amplitude of electric field (the 'wave-breaking limit') that can be sustained in the plasma 
and, as mentioned in the introduction, our aim is to investigate the relationship between 
the shape of the distribution (i.e. E) and the wave-breaking limit. 

For simplicity, we choose to describe the plasma using a distribution / where / = a is 
a positive constant inside a 7-dimensional region hi C £ and / = outside. In particular, 
we consider U to be the union over each point p G of a domain Hp whose boundary 
dWp in £p is topologically equivalent to the 2-sphere. Such distributions are sometimes 
called 'waterbags' in the literature and can be completely characterized by and the 
constant a. Our approach may be considered as a multi-dimensional generalization of the 
purely 1-dimensional relativistic waterbag model used in [5] to examine wave-breaking. 

We work in Minkowski spacetime {A4 , g) and assume that the ions constitute a homo- 
geneous immobile background. We employ an inertial coordinate system (a;") adapted 
to the ions: 

(34) Nion = ?^ionl9o, 

where 

(35) g = -dx° (g) dx° + dx'^ ® dx^ 4- dx^ ® dx'^ + dx^ (g) dx^ 

and the ion proper number density njon is constant. 

To proceed further we seek a form for E axisymmctric about x^ whose pointwise 
dependence in A4 is on the wave's phase C. = x^ — vx^ only, where v is constant and 
< w < 1. We suppose that all electrons described by / arc travelling slower than the 
wave, and the wave 'breaks' if the longitudinal velocity of any plasma electron equals v 
(i.e. an electron 'catches up' with the wave). 

Introduce 

(36) = vdx^ - = dx^ - vdx°, 
and decompose as 

(37) = [m(C) + + V'(C\ C) + i?sin(ei) cos(e')dxi + Rsm(C^) sm{^^)dx^ 
for < ^1 < TT, < < 27r where i? > is constant. 
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Here, (76^, 76^, dx^, c?x^), with 7 = l/yl— t^, is an orthonormal coframe on M 
adapted to C,. Since is future-directed and timelike, and is timclikc, it follows 
ei(Vf) < and/i + ^1(^1) > 0. 

The component ip is determined using (18), 



(38) V = - ^[^l + A{eW-^^[\ + R^sm\e)l 

where the negative square root is chosen because no electron is moving faster along 
than the wave. 

Substituting the ansatz (37) together with a purely longitudinal electric field depend- 
ing only on 

(39) F = E{C)dx" Adx^, 
into (33) yields 



„ 1 m dn 

(40) E^---^. 

7 9 dC 



Equation (40) is used to eliminate E from Maxwell equations (5) and obtain a differential 
equation for fj.. 

The electron number current is calculated using (9): 



(41) Nip)^a[ / 



and (5, 37, 38) yield 



2 J/-9 ~ ^ion7 

7^ dQ'^ m 



(42) 



and 



2 r / \ 1/2 

^2^i?2^ + " "^'[^ + sin2(ei)] j sinie) cos(e') 



(43) 27ri?2 y ^(^1) sin(e^) cos(e^) d( 



1\ Ji-l _ n-.onl'^V 



where a is the value of / inside Wp. 

The form of the 2nd order autonomous non-linear differential equation (42) for /i is 
fixed by specifying the generator A{S}) of dWp subject to the normalization condition 
(43). 
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4'1. Electrostatic wave-breaking. - The form of the integrand m (42) ensures that the 
magnitude of oscillatory solutions to (42) cannot be arbitrarily large. For our model, the 
wave-breaking value /^wb is the largest fi for which the argument of the square root in 
(42) vanishes, 



(44) 



because ^ < ^wb yields an imaginary integrand in (42) for some The positive square 
root in (44) is chosen because, as discussed above, n + A{(^^) > and in particular 
/iwb+A(ei) >0. 

The wave-breaking limit -Emax is obtained by evaluating the first integral of (42) 
between ^wb where E vanishes and the equilibrium(^) value ^eq of /i where ii^ is at a 
maximum. Using (43) to eliminate a it follows that fx^q satisfies 

7r 

i J Aie)smie)cos{e)de 



(45) ^ J (^Wq + A{e)?-l'[l + R'sm^e)]^ sm{e)cos{e)de 



with 



(46) 



TT 

J A(e^)sin(^i)cos(^i)de^ < 



since a,v > 0. Equation (42) yields the maximum value -Bmax of i?. 



Meq TT 



jA{e')smie')(^os{e')de 



1 /2 

I I (^[^l + A{e)f-^■'[l + R^sin\e)?j sm{e)cos{e)ded^^ 



(47) 



Mwb 



The above is a general expression for i?,nax given A{£^^) as data. In the following, we 
determine a simple expression for an upper bound on E-a^g,^ when A[^^) = — acos(^"'^) 
where a is a positive constant (a > ensures (46) is satificd). Using (47) it follows 



(48) £;^ax = 2TOnion < - Mcq + Mwb + ^ ~ / [^+ (m) + 1- ((")] d^ 



(^) Note that the equilibrium of /i need not coincide with the plasma's thermodynamic equilib- 
rium. 
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where 

1 1/2 

(49) X±(m) (^[M±«x]'-7'[l + i?'(l-x')]) Xrfx 



Furthermore, (45) may be written 

(50) ^^[X+(//eq)+I-(Mcq)] =1 

and since X_|-(/Xcq) > I+(/i) and X_(/x„b) > 2^-(m) for Mwb < ^ /^oq, using (48, 50) 

(51) £^^ax < mnion{^^cq " /iwb) [2^- (/iwb) - X-(/ioq)] 



Furthermore — X_(/icq) < ^yMeq ~ 7^ s-^id I_(/Xwb) < so 

(52) i^niax < ^"^(/^-^q ' Mwb)y/ieq " 7^ 



where Wp — yrikm^V^ the plasma angular frequency (in units where Eq = 1 ^md 
c=l). 

4'2. Wave-breaking in an external magnetic field. - In tackling (32), one may opt to 
seek only those satisfying (33); this approach was followed in the preceeding sections. 
Although, at first sight, this method appears to be a simpler than attempting to solve 
(32), it is not always the simplest option. There are potential advantages in considering 
(32) in its generality, as we will now argue. 

The component of the magnetic field parallel to the velocity of a point charge does 
not contribute to the Lorentz force on that point charge. A similar observation may also 
be applied to certain in (32) even though the (i9i, 82) components of are non-zero. 
Furthermore, the results of the previous section are unaffected by a constant magnetic 
field aligned along x^. 

The axially symmetric introduced above is of the general form 



(53) V^ = {l + Y^ + Zy/^do + y cos(e2)9i +Ysmi^^)d2 + Zds, 
where Y = Y{x,£,^) and Z — Z{x,£^^). Suppose F is of the form 

(54) F = Fi + Fn, 

(55) Fi ^ E{C)dx" A dx\ 

(56) Fii = Bdx^ A dx^ 



We have 



(57) 



iv^Fii ^ BY{cos{^^)dx^ - sm{^^)dx^) 
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and furthermore 



(58) — = -Ysm{e), 

f)y'2 

(59) -Q^-Y^^^i^') 
so 

(60) gab^dx^ = -Y sin{^'^)dx^ +Y cos{£,'^)dx'^ . 

ot, 

Hence, using (26) it follows 

(61) —Lv,FiiAn. = 

m ^ 

and from (32) 

(62) (y V ±, Fi)An^^O. 

Therefore, if satisfies (32) with F ~ Fj, the same velocity field also satisfies (32) 
with F = Fj + Fij. It follows that a longitudinal magnetic field does not influence 
an axially symmetric discontinuity in the electron distribution and the results of the 
previous section hold for non-zero constant B. 

5. — Beyond the Lorentz force 

The previous discussion clearly shows that there is merit in considering (32) in its 
generality. We argue that for future extension of this work to fields with more complicated 
spacetime dependence, it is prudent to eschew (33) in favour of (32). We now illustrate 
this point further using a very simple example. 

Let the chain E be such that 

(63) V^ = {1+ R^fl^do + Rs\\i{£})cos{^'^)di + Rsm{S}) s\M{^^)d2 + i? cos(^i)93 
where i? is a function on M and, using (26), 

(64) = i?2 sm{S}) cos{S})dx^ A dx^ + R^ iiui^{S})dx^ A (sin(^2)dxi - cos(C^)di 

The 4-acceleration of is 
^ ~ RV^R , 



(65) + V^r{^ sin(^i) cos(^2)^^i sin(^i) sm{£,'^)dx'^ + coB{^^)dx^ ' 
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and, for simplicity, we assume that (65) is in response to a fongitudinal electric field, 
(66) F^Edx°/\dx^, 
which contributes to (32) as 



(67) Lv.F = £;( (1 + i?2)i/2^a;3 _ Rcos{^^)d. 



Clearly the Lorentz force equation (33) cannot be satisfied for general R, since Vy^V^ 
contains terms in dx^, dx"^ which cannot cancel against terms in iVt^F. However, the only 
nonzero contribution to the left-hand side of (32) can be made to vanish by requiring 



(68) Ffi? = ^Ex/TTb^ cos(^i) 

TO 



Inspection of the ^ dependences in (63) and (68) reveals that R can depend only on 
and 



(69) — v/l + i?2 = ±E 

dx-^ TO 

so E also depends only on x^. 

To compare the above with solutions to (33), we seek a reparametcrisation of S whose 
corresponding family of 4- velocities satisfies (33). In particular, we consider a map p 

p:V xV' xV 

(70) - {^,e = i^'i^^a^e = i>\x,a) 

where V C and (C'\^'^) £ C Then given S satisfying = 0, 

(71) {T.O pfu ^ p*{Y.*uj) ^Q. 

The chains S and (S o p) locally represent the same discontinuity, and are physically 
equivalent. However, the families = V^da and W^i = W^,da of vector fields, where 

(72) = W^, = (S o pfx", 
are different. We demand 

(73) Vw'M' = -iw,,F 



TO 



and using (71) it follows 



(74) Wf'iP- = 0, Wj'?/'^ = --^Vl + sin(V'^). 

mR 

One may solve (74) to determine (Sop), but is clear that (for general E) solving for the 
discontinuity in terms of S is a simpler task. Furthermore, we expect this state of affairs 
to hold for more complicated configurations. 



12 



D.A. BURTON, A. NOBLE, H. WEN 



6. Conclusion 

Wc have developed a covariant formalism for tackling discontinuities in 1-particlc 
distributions. We have used it to develop wave-breaking limits for models of thermal 
plasmas whose distributions have effectively 1-dimensional spacetime dependence but 
are 3-dimensional in velocity. 



Appendix A. 



Vertical and horizontal lifts 

Given tensors on a manifold A4, there are a number of ways to lift them onto the 
tangent manifold TAi. The simplest and best known of these are the vertical lift and 
(given a connection on M) the horizontal lift. The following is a summary of some 
important properties of these lifts; for more details see, for example. [15]. 

The vertical lift is a tensor homomorphism. Acting on forms, it is equivalent to the 
pull-back with respect to the canonical projection map tt : TA4 Ai: 

(A.l) = n*(3, y (3e AM. 

Acting on a vector Y, the vertical lift is 

(A.2) y^ = y°— yY^Y^—eTM. 

Note from these definitions, contraction of the vertical lift of a vector with the vertical 
lift of a form vanishes: 



(A.3) t(yv)/3^ = 0. 



The horizontal lift makes use of a connection V, with coefficients r"5c in a coordinate 
basis da ~ d/dx"^: 

(A.4) Vs^db = T\,da. 

The horizontal lift dx"^^ of the basis form dx" is 
(A.5) dx''" = dx" x^VXJx^ 

while that of the basis vector djdx'^ is 

The horizontal lifts of more general 1-forms and vectors may be determined from the 
relations 

(A.7) (//3)« = /^/3« /?eAiAl, 

(A.8) (/r)^ = jvyjf V r e TTW, 
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for any function / on . 

Similarly to the vertical lift, the contraction of the horizontal lift of a vector with the 
horizontal lift of a form vanishes: 



Two other useful identities relate to the contractions of vertical and horizontal lifts 
of forms and vectors: 
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